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Abstract 

Searches for HZ production with the Higgs boson decaying into an invisible 
final state were performed using the data collected by the DELPHI experiment 
at centre-of-mass energies between 188 GeV and 209 GeV. Both hadronic and 
leptonic final states of the Z boson were analysed. In addition to the search 
for a heavy Higgs boson, a dedicated search for a light Higgs boson down to 
40 GeV/c^ was performed. No signal was found. Assuming the Standard Model 
HZ production cross-section, the mass limit for invisibly decaying Higgs bosons 
is 112.1 GeV/c^ at 95% confidence level. An interpretation in the Minimal 
Supersymmetric extension of the Standard Model (MSSM) and in a Majoron 
model is also given. 
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1 Introduction 

The data collected by DELPHI have been searched for the presence of a Higgs boson 
produced in association with a Z, in e"'"e~ HZ, but which decays to stable non- 
interacting particles. The process is illustrated in Fig. 1. Such invisible Higgs boson 
decays can occur in Supersymmetry, where the Higgs could decay into a pair of neutrali- 
nos Xi [1^3]. In such models Xi is assumed to be the lightest supersymmetric particle 
and therefore assumed to be stable. It is weakly interacting with ordinary matter. Invis- 
ible Higgs decays also occur in Majoron models [4-6] with the Higgs decaying into two 
Majorons. The results of the search described in this article are valid more generally in 
models with stable Higgs bosons that do not interact in the detector. 

Similar searches have been previously performed by DELPHI [7,8] using data at lower 
centre-of-mass energies and by other LEP experiments [9,10]. In this paper searches are 
presented in four different final states, where the Z decays either into a qq, e+e", 
or r+r^ pair. 

e- LSP, J 



LSP, J 



q,l 



e+ q, 1+ 

Figure 1: Feynman diagram describing the HZ production with the Higgs boson decaying 
into invisible particles, e.g. the lightest supersymmetric particle (LSP) or a Majoron (J) 
in models with an extended Higgs sector. 

The paper is organised as follows: First the analyses in the hadronic channel are 
addressed separately in high and low mass ranges. Then we describe the analyses in the 
leptonic channels which cover /i, e, and r final states. Next, the results are summarised 
and 95% Confidence Level (CL) limits are calculated. The limits are then reinterpreted in 
the framework of the Minimal Supersymmetric extension of the Standard Model (MSSM) 
and in a Majoron model. 

2 The DELPHI detector and the data set 

The analyses were mainly based on the information from the tracking system, the 
calorimeters, the muon chambers, and the photon veto counters of the DELPHI detector. 
The scintillation counters veto photons in blind regions of the electromagnetic calorime- 
ters at polar angles near 40°, 90° and 140°. The DELPHI detector and its performance 
are described in detail in Ref. [11,12]. 
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Year 


Low mass range(GeV/c^) 


High mass range(GeV/c^) 


1998 


40-90 


75-120 


1999 


40-100 


75-120 


2000 


40-105 


95-120 



Table 1: Hadronic channel: Low and high Higgs boson mass ranges for three years of 
data-taking. 

The data set analysed in this paper was taken in the years 1998 to 2000. In 1998 and 
1999, data were recorded at centre-of-mass energies 188.7, 191.6, 195.6, 199.6 and 201.7 
GeV. In 2000 the LEP energy was varied from 199.7 to 208.4 GeV and the data taken at 
energies below and above 205.8 GeV were analysed as two independent subsamples, with 
mean energies of 205.0 and 206.7 GeV. At the end of the year 2000 data taking, one of the 
twelve sectors of the Time Projection Chamber (TPC) became non-operational. Data 
taken afterwards were then treated as a separate sample, with a mean centre-of-mass 
energy of 206.3 GeV. In the following, these three subsamples of the 2000 data set will 
be referred to by the energy of each simulation for the corresponding data, namely 205.0, 
206.5 and 206. 5U. The simulation of the last data taking period (206. 5U) included the ef- 
fect of the missing TPC sector in the detector setup and the changes in the reconstruction 
software to partly recover this loss. 

For the analysis of the hadronic and leptonic channels different criteria are required 
on the detector status during data taking. As a result the total data sets correspond to 
589 pb^^and 571 pb~^, respectively. For the simulation of the signal the HZHA generator 
[13] was used for the four final states. For all the years of data-taking simulated signal 
samples with 5000 events per mass point and channel were generated with the Higgs 
masses from 40 to 90 GeV/c^ in 5 GeV/c^ steps, from 90 to 115.0 GeV/c^ in 2.5 GeV/c^ 
steps and at 120 GeV/c^ . 

The background processes e+e~^qq(n7) and fi^ fi~ {n.'-f) were generated using the KK2F 
generator [14] and the background process T^T~{nj) was generated using the KORALZ 
generator [15]. The processes which lead to charged and neutral current four-fermion 
final states were generated with the WPHACT generator [16]. The PYTHIA generator [17] 
was used to describe the hadronic two-photon processes and the BDK generator [18] was 
used to describe the leptonic two-photon processes. Finally, the BHWIDE generator [19] 
was used for the Bhabha processes. Both signal and background events were processed 
through the full DELPHI detector simulation [11]. The inoperative sector in the TPC is 
also taken into account in the corresponding simulation in the 206. 5U data set. 

3 The hadronic channel 

The hadronic decay of the Z represents 70% of the HZ final states. The signature 
of an invisible Higgs boson decay is a pair of acoplanar and acoUinear jets with a di-jet 
mass compatible with the Z mass and missing energy and momentum due to the invisibly 

decaying Higgs boson. 

In order to obtain a good performance in the whole mass range, two overlapping mass 
windows were defined for each year of operation and the analyses were optimised for each 
window as defined in Table 1. 
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3.1 High mass analysis 

The selection of HZ candidate events consists of several steps in order to suppress the 
bulk of the background. First, the events were clustered into jets using the DURHAM [20] 
algorithm. Then a preselection was applied to remove most of the two-photon background 
and a great part of the backgrounds due to four-fermion processes and to hadronic events 
with a radiative return to an on-shell Z. Then the final separation between the signal 
and the background channels was achieved through an Iterative Discriminant Analysis 
(IDA) [21]. The details of the preselection are: 

• Anti-77: Each event was required to have at least 9 charged particle tracks. Two of 
them must have transverse momentum greater than 2 GeV/c and impact parameters 
to the primary vertex less than 1 mm in the transverse plane and less than 3 mm 
along the beam axis. It was also required that the charged energy be greater than 
0.16y^. There should be no electromagnetic shower with more than 0A5^/s, the 
transverse energy^ be greater than 0.15^/s and the sum of the longitudinal momenta 
be greater than 0.25y/s. 

• Anti-qq(n7) and anti-WW: A cut in the ^p^.^ vs. y/s' [22] plane was apphed, required 

40° < . < 140° and 
> 115 GeV 

where y/s' stands for the effective centre-of-mass energy after the initial state radia- 
tion of one or more photons and ^p^^^.^ is the polar angle of the missing momentum. 
In addition, it was required that less than O.OSy^ was deposited in the STIC^ [11], 
v^/x/s was less than 0.96 and that the total electromagnetic energy within 30° of 
the beam directions was less than 0.16i/i. In order to suppress badly reconstructed 
events, candidates in which a jet pointed to the insensitive region between barrel 
and endcap detectors or where both jet axes were below 12° were rejected. A her- 
meticity veto algorithm [23] using the scintillator counters was applied to ensure 
that no photon escaped in the insensitive region of the electromagnetic calorimeter 
at polar angles near 40°, 90° and 140°. To suppress background from WW pair 
production, the energy of the most energetic particle was required to be less than 
0.2y^ and the transverse momentum of any particle in the jet with respect to its 
jet axis (forcing the event into a two-jet configuration) to be less than O.OSy^/c. 
Finally, upon forcing the event into a three-jet configuration, it was required that 
every jet had at least one charged particle in order to suppress qq{n'y) events. 

Twelve variables were used to construct an effective tagging variable in the framework 
of the IDA. In order to calculate these variables, the event was forced into two jets. The 
variables are: 

• E.y/E'^: the normalised energy of a photon assumed to have escaped in the beam 
direction, deduced from the polar angles of the two jet directions in the event. The 
photon energy was normalised to the energy expected for a photon recoiling against 
an on-shell Z. 

• ln(j9T [GeV/c]): logarithm of the transverse momentum of the event. 

• -Evis/v^: visible energy of the event, normalised to the centre-of-mass energy. 

• E^:/ ^/s: transverse energy of the event, normalised to the centre-of-mass energy. 

• ^cone- The minimum polar angle defining a cone in the positive and negative beam 
directions containing 6% of the total visible energy. 

^The transverse energy is the energy perpendicular to the beam axis, defined as = -^pI +Py + 
^ Small angle Tile Calorimeter, covering the very forward region. 
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Variable 


lower cut 


upper cut 






0.90 




i.75 


4.5 




0.15 


0.6 


Pisol/ 


0.008 


0.18 


logiQ( scaled acoplanarity) 


0.3 


2.5 


Thrust 


0.65 


1.0 


In(acollinearity) 


2.0 


4.5 


ln(max(pT )jet [GeV/c]) 


-0.5 


2.50 



Table 2: Tail cuts used in the high mass hadronic analysis. The variables are described 
in detail in section 3.1. 

• I cos 9p^^._^ I : cosine of the polar angle of the missing momentum. 

• £'isoi/ y/s- energy sum between the two cones, defined by half opening angles 5° and 
Q!max around the most isolated particle. The energy sum is then normalised to the 
centre-of-mass energy. The most isolated particle is defined as the particle with 
momentum above 2 GeV/c with the smallest energy sum in the double cone. In the 
momentum interval from 2 to 5 GeV/c, ctmax is set to 60° in order to maximise the 
sensitivity to isolated particles from tau decays in WW —>■ qqVi/ events. An opening 
angle of 25° is used for particles with momenta above 5 GeV/c. 

• Pisoi/v^: momentum of the most isolated particle, as defined above, normalised to 
the centre-of-mass energy. 

• logxo(scaled acoplanarity): The acoplanarity is defined as 180° — A0, where A0 
is the difference in azimuthal angle (in the plane perpendicular to the beam axis) 
between the two jets. In order to compensate for the geometrical instability of the 
acoplanarity for jets at low angles it was multiplied with the angle between the two 
jets. 

• Thrust: thrust value of the event, computed in the rest frame of the visible system. 

• In(acollinearity): logarithm of the acoUinearity (in degrees) of the two-jet system. 

• ln(max(pT [GeV/c])jet): highest transverse momentum of the jet-particles, defined 
by the transverse momentum of any particle in the jet with respect to the jet axis. 

The cuts listed in Table 2 were applied in the tails of the distribution of these variables 
in order to concentrate on the signal region and to avoid long tails in the input variables 
for the IDA. In addition to the cuts listed in Table 2, the number of electrons or muons 
identified by the standard DELPHI algorithms [11] was required to be less than three. 
The agreement between data and simulation is shown in Table 3 and Fig. 2. There is a 
small excess in data over expected background, which is not concentrated in one bin. 

The IDA is a modified Fisher Discriminant Analysis, the two main differences arc the 
introduction of a non-linear discriminant function and iterations in order to enhance the 
separation of signal and background. Two IDA steps were performed, with a cut after 
the first IDA iteration keeping 90 % of the signal efficiency. In order to have two indepen- 
dent samples for the derivation of the IDA function and for the expected performance, 
the signal and background samples were divided in two equally sized samples. As an 
illustration, the distributions of the two IDA variables at ^/s — 206.5 GeV are shown in 
Fig. 3. The shght disagreement in the rates observed at the preselection level is effectively 
removed by the IDA analysis, since it is concentrated mostly outside the signal region. 
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Anti-77 


Anti-qq(n7) & anti-WW 


Tail cuts 


y\jn V ) 


U cl 1.1 ft 










ATP 


ioo.U 


15115 


1 4967 + 81 


1 578 


1 565 2 + 6 


494 


485 9 + 3 2 


1 Ql 


2394 


2351.8 ± 1.3 


258 


249.9 ± 0.9 


88 


79.0 + 0.5 


195.5 


7040 


6782.4 ±3.7 


739 


734.9 ±2.8 


242 


242.0 ± 1.4 


199.5 


7296 


7168.9 ± 3.9 


781 


795. 1 ± 2.8 


295 


261.1+ 1.6 


201.6 


3557 


3407.8 ± 1.9 


396 


382.9 ± 1.3 


152 


130.6 ±0.7 


205.0 


6272 


6011.6 ±3.7 


678 


686.2 ±2.4 


240 


239.2 ± 1.3 


206.5 


6772 


6697.0 ±4.5 


798 


768.5 ± 2.9 


283 


268.2 ± 1.6 


206.5U 


4472 


4560.4 ± 3.9 


534 


541.5 ±2.6 


202 


190.7 ± 1.5 



Table 3: Comparison of simulation and data after the different steps of the preselection 
in the high mass hadronic analysis. The listed errors are from Monte Carlo statistics 
only. 

The observed and expected rates at \fs = 206.5 GeV are shown in Fig. 4 as a function 
of the efficiency to detect a 105 GeV/c^ Higgs boson when varying the cut on the second 
IDA variable. The final cut on the second IDA variable was determined by maximising 
the expected exclusion power. This was done separately for each centrc-of-mass energy to 
optimise the analysis for a 85 GeV/c^ Higgs boson at 188.6 GcV, for a 95 GeV/c^ Higgs 
boson at 191.6 and 195.6 GeV, for a 100 GeV/c^ Higgs at 199.5 and 201.6 GeV and for a 
105 GeV/c^ Higgs at 205.0, 206.5 GeV and 206.5U GeV. Here we assume SM production 
cross-section and a branching ratio of 100% into invisible final states. For example, in 
Fig. 4 the cut on the second IDA is indicated by the dashed vertical line. The final 
number of selected events in data and Monte Carlo simulations is given in Table 7. 

3.2 Low- mass analysis 

For the low-mass analysis, the preselection was adapted for the different event shape 
and kinematics. In the anti-qq(n7) and anti-WW selection the cut in the 9p^.^ vs y/s" 
plane and the cut on \/s'/ y/s were removed in order to increase the signal efficiency. This 
was possible because the signal events have a much smaller amount of missing energy than 
the events in the high-mass range. Some tail cuts were also slightly changed as shown in 
Table 4 and a cut requiring the visible mass to be at least 20% of -^/s was added. Figure 5 
and Table 5 show the agreement of data and background at the preselection level. Figure 
5 a) shows an excess of data over the expected background near due to an 

underestimation of the two-fermion processes. This region is effectively removed by the 
IDA analysis. 

The low-mass analysis also used two IDA steps in order to obtain optimal signal to 
background discrimination. The distributions of the two IDA variables at v^=195.5 GeV 
are shown in Fig. 6. The observed and expected rates at ^/s — 195.5 GeV are shown 
in Fig. 7 as a function of the efficiency to detect a Higgs boson when varying the cut 
on the second IDA variable. The cut on the second IDA variable was again determined 
separately for each centre-of-mass energy as described above. It was optimised for a 
60 GeV/c^ Higgs boson mass at all energies. The final number of selected events in data 
and Monte Carlo simulations is given in Table 8. 
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Variable 


lower cut 


upper cut 






1.20 






0.6 


Pisoi/ y/s 




0.18 


logio (scaled acoplanarity) 


1.0 


2.5 


In(acollinearity) 


2.25 


4.5 



Table 4: Tail cuts used in the low mass hadronic analysis. The variables are described in 
detail in section 3.1. 



(GeV) 


Data 


Anti-77 
MC 


Anti-q 
Data 


q(n7) & anti-WW 
MC 


T 

Data 


ail cuts 
MC 


188.6 


15115 


14967.0 ±8.1 


6604 


6735.2 ± 11.0 


622 


652.0 ±3.9 


191.6 


2394 


2351.8 ± 1.3 


1013 


1051.2 ±1.7 


112 


103.0 ±0.6 


195.5 


7040 


6782.4 ±3.7 


2939 


3003.0 ±4.8 


322 


301.3 ± 1.8 


199.5 


7296 


7168.9 ±3.9 


3122 


3117.7 ±5.0 


338 


315.1 ± 1.8 


201.6 


3557 


3407.8 ± 1.9 


1551 


1495.9 ±2.4 


168 


152.1 ±0.8 


205.0 


6272 


6011.6 ±3.7 


2617 


2614.9 ±4.3 


344 


307.3 ± 1.6 


206.5 


6772 


6697.0 ±4.5 


2885 


2909.0 ±5.2 


305 


293.6 ± 1.7 


206.5U 


4472 


4560.4 ± 3.9 


1878 


1982.5 ±4.7 


257 


237.5 ± 1.7 



Table 5: Comparison of simulation and data after the different steps of the preselection 
in the low mass hadronic analysis. The errors given are from Monte Carlo statistics only. 

3.3 Mass reconstruction 

The recoil mass to the di-jet system corresponds to the mass of the invisible Higgs 
boson. It was calculated with a Z mass constraint for the measured di-jet system from 
the visible energy Eyis and the visible mass myis- The following expression was used 

where pmis is the missing momentum and mz is the Z mass. The recoil mass distribution 
after the final selection for the high-mass analysis is shown in Fig. 8. For the low-mass 

region this method was also used. In cases where the fit obtained negative mass squares 
the standard missing mass calculation a/ -E^is — p^is '^^s used, where i?mis = -s/s — Eyig. 
The recoil mass distribution for the low mass analysis is shown in Fig. 9. 

3.4 Systematic errors 

Several sources of systematics have been considered, first the effect of modelling the 
qq(n7) background from different generators was studied by replacing the KK2F generator 
with the ARIADNE generator [24] at 206.5 GeV. The results were identical within statistical 
errors. The error of the luminosity is conservatively taken to be ±0.5%. The process 
e+e" — qqez/ provides about a fifth of the background and the uncertainty on the cross- 
section of this process is taken to be ±5% [25]. This leads to an ±1% uncertainty of the 
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background. In order to see the influence of tlie jet clustering algorithm the DURHAM 
algorithm was replaced by the LUCLUS algorithm [26]. This results in an uncertainty 
on the background estimation and the signal efficiency in the order of ±1% for the high 
mass regime and an error in the order of ±2.5% for the low mass regime. 

The data and Monte Carlo simulation were found to be in good overall agreement. 
However, since we are searching for events with a large amount of missing energy, we 
become sensitive to the tails of the distributions from the expected Standard Model 
background events. When analysing the same topology for the measurement of the Z 
pair production cross-section [27], it was found that the small disagreement in the tails 
can be cured if the particle multiplicities of data and Monte Carlo simulation are brought 
into agreement. In order not to bias the present analysis, where the disagreements in 
the tails could come from new physics, the tuning of the particle multiplicities was done 
with Z — >■ qq events taken at \/s = 91.1 GeV for each year of data taking. The particle 
multiplicities were estimated separately for the barrel (cos^ < 0.7) and the forward region 
(cos^ > 0.7) and for different momentum bins and separately for neutral and charged 
particles. For each of these classes of multiphcities a separate correction factor P was 
calculated using 

p^ < iVdata > - < NmC > 

< Nmc > 

where < A^^ata > is the mean value of the particle multiplicity in the data and < Nmc > 
is the corresponding simulated value. These correction factors are of the order of a few 
percent in the barrel region, they tend to be larger in the forward region and are also 
larger for neutral than for charged particles. The correction factors obtained were then 
applied to the high energy LEP2 Monte Carlo simulation on an event by event basis. 
The factor P was used as a probabihty to modify the particle multiplicities in the Monte 
Carlo simulation and related variables were recalculated. If P was less than zero, there 
were fewer particles in data than in Monte Carlo and the particles of the corresponding 
class were removed in the simulated events. For P greater than zero, particles have to be 
added to the simulated events. This was performed copying another particle of the same 
class and smearing its momentum by 2.5% in order not to affect the event jet topology. 
If there was no particle of the corresponding class, a particle of the adjacent class was 
taken and scaled to fit into this class. Note that these modifications of the multiplicities 
in the Monte Carlo simulation were not used to change the analysis, but only to estimate 
the systematic errors. The effect on the final background estimation ranges from ±10.5% 
(1998), ±4.7% (1999) to ±10.6%(2000) for the high mass range analysis. For the low 
mass range the effects are smaller, they range from ±6.6 % (1998), ±4.3% (1999) to 
±5.6% (2000). This procedure also affects the signal efficiencies leading to a reduction 
of the relative signal efficiency of up to ±1.5%. The application of this method to the 
analysis variables leads to a better agreement of data at the preselection level as has 
been observed previously in the measurement of the Z pair production cross-section [27], 
leading to a better estimation of the systematic error on the simulated background. The 
total systematic error and statistical error from the limited MC statistics are combined 
in quadrature and given in Table 7 and Table 8. 

4 Leptonic channels 

The leptonic channel represents about 10% of the HZ final state. The experi- 

mental signature of the HZ(Z £^i^) final states is a pair of acoplanar and acollinear 
leptons, with an invariant mass compatible with that of an on-shell Z boson. 
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The analysis contains a preselection for leptonic events. Then, the search channel 
is defined by the lepton-type of the Z decay mode and for each decay mode specific 
selection cuts were applied. Two different sets of final cuts were used, depending on the 
reconstructed mass, defining the low-mass and high-mass ranges. 

4.1 Leptonic preselection 

To ensure a good detector performance the data corresponding to runs in which sub- 
dctcctors were not fully operational were discarded. In particular it was required that the 
tracking subdetectors and calorimeters were fully operational and that the muon cham- 
bers were fully functional. This resulted in slightly smaller integrated luminosities than 
for the hadronic search channel. An initial set of cuts was applied to select a sample 
enriched in leptonic events. A total charged-particlc multiplicity between 2 and 5 was 
required. All particles in the event were clustered into jets using the LUCLUS algo- 
rithm [26] ((ijoin = 6.5 GeV/c) and only events with two reconstructed jets were retained. 
Both jets had to contain at least one charged particle and at least one jet had to contain 
exactly one charged particle. 

In order to reduce the background from two-photon collisions and radiative di-lepton 
events, the acoplanarity, 6*300?, had to be larger than 2°, and the acoUinearity, Og^oi, had to 
be larger than 3°. In addition, the total momentum transverse to the beam direction, px, 
had to exceed 0.02y^/c. Finally, the energy of the most energetic photon was required to 
be less than 0.15y/s. The angle between that photon and the charged system projected 
onto the plane perpendicular to the beam axis had to be less than 170°. The agreement 
of data and background at the preselection level is shown in Fig. 10 for all data sets. 



(GeV) 


Data 


UC 


Data 


e+e 

MC 


Data 


MC 


188.6 


64 


49.7 ±0.8 


314 


298.0 ±7.1 


124 


148.2 ±3.7 


191.6 


10 


7.9 ±0.1 


46 


45.5 ± 1.1 


18 


22.4 ±0.8 


195.5 


19 


22.5 ±0.3 


132 


125.5 ±3.1 


78 


62.2 ±2.1 


199.5 


24 


24.8 ±0.3 


149 


134.6 ±3.2 


81 


74.9 ± 1.8 


201.6 


17 


12.1 ± 0.2 


60 


65.9 ± 1.6 


31 


32.1 ± 1.1 


205.0 


11 


20.5 ±0.3 


98 


114.2 ±2.7 


70 


55.3 ± 1.0 


206.5 


26 


23.1 ±0.3 


110 


129.5 ±2.2 


76 


71.9 ± 1.1 


206.5U 


6 


14.6 ±0.2 


79 


76.3 ± 1.9 


48 


40.3 ± 1.3 



Table 6: Comparison of simulation and data at preselection level in the three leptonic 
channels. The errors refiect the Monte Carlo statistics only. The last fine (206. 5U) refers 
to the data taken with one TPC sector inoperative, which has been fully taken into 
account in the event simulations. 



4.2 Channel identification 

For the preselected events, jets were then identified as either /x, e or r and two leptons 
with the same fiavour were required. Owing to the low level of background, the three 
lepton identifications rely on loose criteria. A charged particle was identified as a muon if 
at least one hit in the muon chambers was associated to it, or if it had energy deposited 
in the outermost layer of the hadron calorimeter. In addition, the energy deposited in 
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the other layers had to be compatible with that from a minimum ionising particle. Only 
jets with exactly one charged particle were tagged as muons. For the identification of a 
charged particle as an electron the energies deposited in the electromagnetic calorimeters, 
in the different layers of the hadron calorimeter, and in addition the energy loss in the 
Time Projection Chamber were used. An electron jet had to contain a maximum of 
two charged particles with at least one identified electron. A lepton was defined as a 
cascade decay coming from a r if the momentum was lower than 0.13y^/c. In this case 
the charged particle is no longer classified as a muon or as an electron. If no muon or 
electron was identified, the particle was considered a hadron from a r decay. Thus, there 
is no overlap between the event samples selected in the three channels. The number of 
data and simulated background events are given in Table 6 for each centre-of-mass energy. 
A detailed description of the lepton identification is given in Ref . [28] . 

4.3 Channel-dependent criteria 

After the preselection, different cuts were applied in each channel in order to reduce the 
remaining background. The optimisation of the efficiency has been performed separately 
for mass ranges of 50 to 85 GeV/c^ and 85 to 115 GeV/c^. 

In the channel only events with exactly two charged particle tracks were ac- 

cepted. The direction of the missing momentum had to deviate from the beam axis by 
more that 18° in order to reject Z — > (7) and 77 — > processes. The di-muon 

mass was required to be between 75 GcV/c^ and 97.5 GcV/c^, to be consistent with the 
Z boson mass. After that, two different sets of cuts were applied depending on the recon- 
structed Higgs boson mass as defined in section 4.5. If the reconstructed mass was higher 
than 85 GeV/c^ the momentum of the most energetic muon had to be between 0.2i/i/c 
and 0.4y^/c. Furthermore, E^ig < 0.55-v^, px < 0.25^/s/c and ^^acoi < 60° was required. 
Otherwise, the momentum of the most energetic muon had to be between 0.25y^/c and 
0.45v^/c, and 0.45v^ < E'vis < 0.65^, px < 0.4v^/c and 45° < 6'acoi < 85° was 
required. The mass resolution for Z is about 4.5 GeV/c^. 

In the e"'"e~ channel a maximum of four tracks were required. The most important 
background arises from radiative Bhabha scattering and Ze"'"e~ events. To suppress these 
backgrounds, the direction of the missing momentum and the polar angle of both leptons 
had to deviate from the beam axis by more than 18°, the transverse energy had to be 
greater than 0.15i/i and the neutral electromagnetic energy had to be less than 0.1i/i. 
The invariant mass of the two leptons had to be between 75 GeV/c^ and 100 GeV/c^to 
be consistent with the Z boson mass. The mass resolution for Z e+e" is about 
5.7 GeV/c^. Then, if the mass reconstructed was higher than 85 GeV/c^, the momentum 
of the most energetic electron had to be lower than 0.35y^, and the total associated 
energy was required to be less than 0.55-y/s, px < 0.25-v/s/c and ^acoi < 60°. Otherwise, 
the momentum of the most energetic electron had to be between 0.25y^/c and 0.45y^/c, 
and the total associated energy was required to be less than 0.65y^. In addition, the 
selection px < 0.4y^/c and 45° < 6'acoi < 85° was applied. 

In the r+T~ channel tighter cuts were applied on the acoplanarity and acoUinearity in 
order to reduce the remaining backgrounds from r+r~(7) and 77 ii processes. The 
invariant mass of both jets had to be less than 3 GcV/c^. In addition, the transverse 
energy had to be greater than O.ly^, the visible energy of all particles with | cos 6' |< 0.9 
had to be greater than 0.06i/i and the energy of both jets had to be less than 0.26^/s. 
Finally, if the mass reconstructed was higher than 85 GcV/c^, the acoUinearity had to 
be between 10° and 60°, otherwise, it had to be between 45° and 85°. No cut on the 
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reconstructed mass is applied because of the large missing energy from the associated 
neutrinos. 

4.4 Systematic uncertainties 

Several sources of systematic uncertainties were investigated for their effect on the sig- 
nal efficiency and the background rate. The particle identification method was checked 
with di-lepton samples both at Z peak and high energy, and the simulation and data 
rates were found to agree within ±1%. The modelling of the preselection variables agrees 
within statistical errors with the data. The track selection and the track reconstruction 
efficiency were also taken into account in the total systematic error. The effects of detec- 
tor miscalibration and deficiencies were investigated using //+//~7 or e"^e~7 events, where 
the lepton energies arc determined directly and recoiling from the photon. The compar- 
ison between data and simulation rate was found to be better than ±1%. Additional 
systematic effects were estimated by comparing the data collected, at the Z peak, during 
the period with one TPC sector inoperative with simulation samples produced with the 
same detectors conditions. The total systematic error on the signal efficiency was ±1.1%. 
The total systematic error on the background rate was up to 10%. The total systematic 
error and statistical error from the limited MC statistics are combined in quadrature and 
given in Table 7. 

4.5 Mass reconstruction 

The mass of the invisibly decaying particle was computed from the measured energies 
assuming momentum and energy conservation. To improve the resolution a fit was 
applied constraining the visible mass to be compatible with a Z. In the case of the t+t~ 
channel, the measured four-momenta of the decay products do not reproduce correctly 
the r energy. Therefore, the mass was calculated under the assumption that both r 
leptons had the same energy and the r neutrino went along the direction of the r lepton. 
This, together with the visible mass constraint, allowed an estimation of the r energy and 
of the invisible mass. The invisible mass for the candidates as well as for the expected 
background from Standard Model processes for the different channels is shown in Fig. 11. 

5 Results 

A comparison of the observed and predicted numbers of selected events for the four 
channels is summarised in Tables 7 and 8. The agreement between the data and the 

SM prediction is good for all channels and no indication for a Higgs boson decaying 
into invisible particles has been observed. The signal efficiencies of the four channels are 
shown in Fig. 12 as a function of the Higgs mass for = 206.5 GeV. 

5.1 Model independent limits 

The cross-section and mass limits were computed at the 95% CL with a likelihood 
method [29]. One-dimensional distributions of the reconstructed mass serve as input for 
the likelihood calculation. The impact of the correlation of the systematic errors is small 
and the limits result largely from the data taken at the higher centre-of-mass energies. 
More details about the confidence definition and computation can be found in Ref. [23]. 
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Channel 


Luminosity 


Data 


Expected 


Signal efficiency 


(GeV) 








background 


(%) 


188.6 


qq 


152.4 


65 


71.3 ±7.7 


40.9 ± 1.9 


191.1 


qq 


24.7 


2 


5.6 ±0.3 


39.6 ± 1.7 


195.5 


qq 


74.3 


21 


18.7 ± 1.0 


50.8 ± 1.7 


199.5 


qq 


82.2 


21 


20.1 ± 1.0 


51.9 ± 1.7 


201.6 


qq 


40.0 


11 


10.8 ±0.5 


50.7 ± 1.7 


205.0 


qq 


74.3 


9 


12.2 ± 1.3 


36.4 ±2.1 


206.5 


qq 


82.8 


13 


13.5 ± 1.5 


37.0 ±2.1 


206.5U 


qq 


58.0 


11 


8.4 ±0.9 


31.6 ±2.1 


188.6 


11+11- 


153.8 


7 


6.9 ±0.6 


44.0 ± 1.9 


191.1 




24.5 


4 


1.1 ±0.1 


52.8 ± 1.6 


195.5 


fj,+fj,~ 


72.4 


3 


3.5 ±0.2 


63.8 ± 1.5 


199.5 




81.8 





3.9 ±0.3 


63.0 ± 1.5 


201.6 




39.4 


2 


1.8 ±0.2 


62.5 ± 1.5 


205.0 




69.1 





3.0 ±0.3 


62.8 ± 1.5 


206.5 




79.8 


2 


3.3 ±0.3 


62.1 ± 1.5 


206.5U 




50.0 





2.2 ±0.2 


56.9 ± 1.6 


188.6 


e+e" 


153.8 


4 


7.9 ±0.7 


34.2 ± 1.3 


191.1 


e+e" 


24.5 


1 


1.2 ±0.2 


40.8 ± 1.6 


195.5 


e+e" 


72.4 


4 


4.7 ±0.5 


45.3 ± 1.6 


199.5 


c+c^ 


81.8 


5 


4.1 ±0.4 


45.2 ± 1.6 


201.6 


e+e^ 


39.4 


1 


1.9 ±0.2 


45.1 ± 1.6 


205.0 


e+e~ 


69.1 


3 


3.6 ±0.3 


44.8 ± 1.6 


206.5 


e+e" 


79.8 


1 


4.0 ±0.4 


42.9 ± 1.6 


206. 5U 


e+e" 


50.0 


1 


2.3 ±0.3 


39.9 ± 1.6 


188.6 


T+T~ 


153.8 


7 


9.4 ±0.8 


21.4 ± 1.4 


191.1 


T+T^ 


24.5 


1 


1.9 ±0.2 


17.3 ± 1.4 


195.5 


T+T~ 


72.4 


7 


5.7 ±0.6 


20.2 ±2.1 


199.5 


T+T~ 


81.8 


10 


6.3 ±0.6 


27.3 ± 1.5 


201.6 


T+T~ 


39.4 


2 


3.3 ±0.4 


28.2 ± 1.5 


205.0 


T+T~ 


69.1 


5 


5.7 ±0.6 


29.5 ± 1.5 


206.5 


T+T" 


79.8 


3 


7.1 ±0.7 


30.3 ± 1.5 


206.5U 


T+T~ 


50.0 


2 


4.5 ±0.4 


29.5 ± 1.5 



Table 7: Integrated luminosity, observed number of events, expected number of back- 
ground events and signal efficiency (100 GcV/c^ signal mass) for different energies. The 
last lines of each channel (206. 5U) refers to the data taken with one TPC sector inoper- 
ative, which has been fully taken into account in the event simulations. Systematic and 
statistical errors are combined in quadrature from the results of each analysis. 
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1 — 


Channel 


Luminosity 


Data 


Expected 


Signal efficiency 


(GcV) 




/I 1 \ 
(pb ^) 




Background 


(%) 


188.6 


qq 


152.4 


58 


51.5 ± 3.8 


49.1 ± 1.6 


191.0 


qq 


2A.7 


6 


10.1 ± 0.5 


50.0 ± 1.7 


195.5 


qq 


74.3 


36 


31.3 ± 1.6 


49.6 ± 1.7 


199.5 


qq 


82.2 


37 


44.3 ±2.3 


50.5 ± 1.7 


201.6 


qq 


40.0 


10 


12.0 ±0.6 


44.2 ± 1.7 


205.0 


qq 


74.3 


26 


26.2 ± 1.7 


47.0 ± 1.5 


206.5 


qq 


82.8 


30 


33.4 ±2.1 


48.8 ± 1.5 


206.5U 


qq 


58.0 


10 


18.0 ± 1.2 


43.6 ± 1.5 



Table 8: Integrated luminosity, observed number of events, expected number of back- 
ground events and signal efficiency (60 GeV/c^ signal mass) for different energies in the 
low mass analysis. The last lines of each channel (206. 5U) refers to the data taken 
with one TPC sector inoperative, which has been fully taken into account in the event 
simulations. Systematic and statistical errors are combined in quadrature. 

All search channels and centre-of-mass energies were treated as separate experiments to 
obtain a likelihood function. In total 40 channels were evaluated as listed in Tables 7 
and 8, in addition to the qq channels from 161 and 172 GeV data [7], and the qq and 
channels from 183 GeV data [8]. In order to address the overlap between the low and 
high mass analyses in the hadronic channel, the expected performance was calculated for 
both analyses in the overlap region. At each test mass the analysis with the best expected 
exclusion power was then chosen for the calculation of the limit. 

No indication of a signal is observed above the background expectation. This is shown 
in Fig. 13 which displays the curves of the confidence levels in the background hypothesis, 
CLb, as a function of the Higgs boson mass hypothesis, for each channel separately. Over 
most of the range of masses the agreement between data and the background expectations 
is within one standard deviation. However, at a few masses in the muon and electron 
channels, there are disagreements near or slightly above two standard deviations, which 
are due to deficits of data in several bins of the reconstructed mass spectra in these 
channels, as shown in Fig. 11. Figure 14 displays the observed and expected upper limits 
on the cross-section for the process e+e^ — > Z (anything) H (invisible) as a function of the 
Higgs boson mass. From the comparison with the Standard Model (SM) Higgs boson 
cross-section the observed (expected median) mass limits are 112.1 (110.5) GeV/c^ for 
the Higgs boson decaying into invisible particles. 

In a model-independent approach the branching ratio into invisible particles -Bi?inv 
can be considered a free parameter. The remaining decay modes are then visible and 
are assumed to follow the SM decay probabihties. In this case the searches for visible 
and invisible Higgs boson decays can be combined to determine the excluded region in 
the -B-Rinv versus mn plane assuming SM production cross-sections. Using the DELPHI 
data from the SM Higgs searches [7,23,30-32] a lower mass limit of 111.8 GeV/c^ can be 
set independently of the hypothesis on the fraction of invisible decay modes, as shown 
in Fig. 15. In computing these limits, the overlap between the standard Huu and the 
invisible Higgs boson hadronic selections have been avoided, conservatively for the limit, 
by omitting the Rui? (Hjnyqq) results in the region BRi,^y > 50%(< 50%). 
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5.2 Limits for a Majoron model 

The limits computed above can be used to set a limit on the Higgs bosons in a Majoron 
model [4-6] with one complex doublet and one complex singlet rj. Mixing of the real 
parts of (j) and rj leads to two massive Higgs bosons: 

H = (pR cos 6 — TjR sin 6 

S = (pRsinO + Tjji cos 6 

where is the mixing angle. The imaginary part of the singlet is identified as the Majoron. 
The Majoron is decoupled from the fermions and gauge bosons, but might have a large 
couphng to the Higgs bosons. In this model the free parameters are the masses of H and 
S, the mixing angle 6 and the ratio of the vacuum expectation values of the two fields 

and 7] (tan/? = The production rates of the H and S are reduced with respect to the 

SM Higgs boson, by factors of cos^ 6 and sin^ 6, respectively. The decay widths of the H 
and S into the heaviest possible fermion-antifermion pair are reduced by the same factor 
and their decay widths into a Majoron pair are proportional to the complementary factors 
{cos^d for S and sin^ ^ for H). The HZ and SZ cross-section times branching ratio into 
invisible decays is calculated and compared to the excluded cross-section of section 5.1. 
In the case where the invisible Higgs boson decay mode is dominant (tan/? larger than 
about 10), the excluded region in the mixing angle versus Higgs boson mass plane is 
shown in Fig. 16. 

5.3 Limits in the MSSM 

In the MSSM, there are parameter regions where the Higgs boson can decay into 
neutralinos, which leads to invisible Higgs decays. As an illustration a benchmark 
scenario including such decays was defined from the so-called "m^-max scenario" [23] . In 
this scenario the MSSM parameters arc the mass of the pseudoscalar Higgs boson, mA, 
the ratio of the vacuum expectation values, tan/?, the mixing in the scalar top sector Xt, 
the gaugino mass M2 and the Higgs self-coupling ji. M2 and // were modified to obtain 
light neutralino masses setting M2 = n = 150 GeV/c^. Then, a scan was performed in 
the tan l3-mp^ plane. For each scan point the hZ production cross-section and the Higgs 
boson branching ratio into neutralinos were calculated, and the point was considered as 
excluded if the product was found to be larger than the excluded cross-section as shown 
in Fig. 14. Figure 17 shows the excluded region from the search for invisible Higgs decays, 
the theoretically forbidden region, and the region where the branching ratio h — * X^x" is 
less than 1%. In this benchmark scenario, the invisible Higgs boson search covers a large 
region in the low tan (3 regime. The white regions cannot be excluded by the invisible 
Higgs searches alone because the branching ratio into neutralinos is too small. The search 
for the invisible Higgs boson decays also sets limits in the general framework searches 
for Supersymmetric particles [33] and for searches in Anomaly Mediated Supersymmetry 
Breaking (AMSB) models [34]. 

6 Conclusion 

In the data samples collected by the DELPHI detector at centre-of-mass energies from 
189 to 209 GeV, 153 qq (213 for the low mass analyses), 18 jJi^jJT, 20 e+e" and 37 t+t~ 
events were selected in searches for a Higgs boson decaying into invisible final states. 
These numbers are consistent with the expectation from SM background processes. 
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We set a 95% CL lower mass limit of 112.1 GeV/c^ for Higgs bosons with a Standard 
Model cross-section and with 100% branching fraction into invisible decays. Excluded 
parameter regions are given in a simple Majoron model. The invisible Higgs boson 
search is important to cover some parameter regions in the MSSM where Higgs decays 
into neutralinos are kinematically allowed. 
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Figure 2: Hadronic channel high mass analysis: Distribution of the four IDA input 
variables after the final preselection as described in section 3.1: a) E^/E"^; b) ln{p-x) 
in GeV/c; c) In(acollinearity); d) Thrust. 
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Figure 3: Hadronic channel high mass analysis: Distributions for the IDA variables after 
first (a) and second IDA step (b) at y/s = 206.5 GeV. The dashed line indicates the cut 
on the IDA variable. The white histogram shows the expectation of a 105 GeV/c^ Higgs 
signal where the signal rate is enhanced by a factor 20 for (a) and 4 for (b). 
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Figure 4: Hadronic channel high mass analysis: Data and expected background for the 
206.5 GeV centre- of- mass energy as a function of the efficiency for an invisibly decaying 
Higgs boson of 105 GeV/c^. The lines represent the most important backgrounds with 
the sohd black line showing the sum of all the background processes. In addition the grey 
line shows the expectation for a 105 GeV/c^ Higgs signal added on top of the background. 
The vertical dashed line indicates the final cut chosen to maximise the sensitivity. 
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Figure 5: Hadronic channel low mass analysis: Distribution of the four IDA input vari- 
ables after the final preselection as described in section 3.1: a) E^/E'^; b) In(pT) in GeV/c; 
c) In(acollinearity); d) Thrust. 
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Figure 6: Hadronic channel low mass analysis: Distributions for the IDA variables after 
first (a) and second IDA step (b) at a/s=195.5 GeV. The dashed line indicates the cut 
on the IDA variable. The white histogram shows the expectation of a 60 GeV/c^ Higgs 
signal where the signal rate is enhanced by a factor 5 for (a). 
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Figure 7: Hadronic channel low mass analysis: Data and expected background for the 
195.5 GeV centre-of-mass energy as a function of the efficiency for an invisibly decaying 
Higgs boson of 60 GeV/c^. The lines show number of events from the most important 
background reactions and the solid black line shows the sum of all the background pro- 
cesses. In addition the grey line shows the expectation for a 60 GeV/c^ Higgs signal 
added on top of the background. The vertical dashed line indicates the final cut chosen 
to maximise the sensitivity. 
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Figure 8: Hadronic channel high mass analysis: Reconstructed Higgs boson mass for 
^/s from 189 to 209 GeV after the final selection. The white histogram corresponds to 
a Higgs boson with 100 GeV/c^ mass decaying with a branching fraction of 100% into 
invisible modes. 
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Figure 9: Hadronic channel low mass analysis: Reconstructed Higgs boson mass for ^/s 
from 189 to 209 GeV after the final selection. The white histogram corresponds to a Higgs 
boson with 60 GeV/c^ mass decaying with a branching fraction of 100% into invisible 
modes. 
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Figure 10: Leptonic channels: Acollinearity distribution for -^i from 189 to 209 GeV 
after the preselection. 
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Figure 11: Leptonic channels: Reconstructed Higgs boson mass in (a) the HyU+yU~channel, 
(b) the He"*" e~ channel and (c) the Hr+r~channel for 189 to 209 GeV after the final 
selection. The white histogram corresponds to a Higgs boson with 100 GeV/c^ mass 
decaying to 100% into invisible modes. 
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Figure 12: Efficiencies for the Higgs boson masses between 40 and 120 GeV/c^ for the 
different selection channels at ^/s = 206.5 GeV. 
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Figure 13: Confidence levels for the different decay channels as a function of the Higgs 
mass. Shown are the observed (solid) and expected (dashed) confidences for the back- 
ground-only hypothesis in the Hqq(a), H/i"'"/i~(b), He"'"e~(c) and Hr"'"r~(d) channels. The 
dark grey band corresponds to the 68.3% expected confidence interval and the light grey 
band to the 95.0% confidence interval. The structures near 94 and 96 GeV in plot (a) are 
due to the switching from the low-mass to the high-mass optimization in the hadronic 
channel. 
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Figure 14: The 95% CL upper limit on the cross-section e"'"e~ — > Z(anything) H(invisible) 
as a function of the Higgs boson mass. The dashed line shows the standard model 
cross-section for the Higgs boson production with BR^^^ — 1. 
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Figure 15: The Higgs boson mass limits as a function of the branching ratio into invisible 
decays Si?inv, assuming a 1 — BRi^y branching ratio into standard visible decay modes. 
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Figure 16: Limit on sin^ 6' as a function of the Higgs boson mass at 95% CL. S and H 
are the Higgs bosons in the Majoron model. The grey region is excluded for the S Higgs 
boson and the hatched region for the H Higgs boson. The massive Higgs bosons decay 
almost entirely into invisible Majoron pairs for large tan/3 values. 
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Figure 17: Excluded region in the MSSM from searches for a Higgs boson decaying into 
invisible final states for the modified "mh-max scenario" described in the text. The 
different grey areas show the theoretically forbidden region (dark), the region where the 
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